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1. Introduction 

Nonlinear differential-difference equations are always more important in applica- 
tions. They enter as models for many biological chains, are encountered frequently in 
queuing-problems and as discretisations of field theories. So, both as themselves and 
as approximations of continuous problems, they play a very important role in many 
fields of mathematics, physics, biology and engineering. Not many tools are available 
to solve such kind of problems. Apart for a few exceptional cases the solution of non- 
linear differential-difference equations can be obtained only by numerical calculations or 
by going to the continuous limit when the lattice spacing vanishes and the system is 
approximated by a continuous nonlinear partial differential equation. Exceptional cases 
are those equations which, in a way or another, are either linearisable or integrable via 
the solution of an associated spectral problem on the lattice. In all such cases we can 
write down a denumerable set of exact solutions corresponding to symmetries of the non- 
linear differential-difference equations. Such symmetries can be either depending just on 
the dependent field and independent variable and are denoted as point symmetries or 
can depend on the dependent field in various positions of the lattice and in this case 
we speak of generalised symmetries. Any differential-difference equation can have point 
symmetries, but the existence of generalised symmetries is usually associated only to the 
integrable ones. 

Few classes of integrable nonlinear differential-difference equations are known [1-4] 
and are important for all kind of applications both as themselves and as starting point 
for perturbation analysis. However not all cases of physical interest are covered and so 
it would be nice to make the classification of integrable nonlinear differential-difference 
equations to provide a basic list of equations to be used as model of nonlinear systems on 
the lattice or as starting point of perturbation theory. A way to accomplish such a goal 
is provided by the so-called formal symmetry approach introduced by A.B. Shabat and 
collaborators in Ufa (see e.g. review articles [5-7]) which consists in classifying all equa- 
tions of a certain class which possess few generalised symmetries. Such an approach has 



1 



been introduced at first to classify partial differential er^uations but then the procedure 
has been extended to the case of differential-difference equations [2-4]. 

In such an approach, one provides conditions under which one can prove the existence 
of at least one (or more) generalised symmetries. These conditions arc basic tools to start 
the procedure of classification, i.e. we look for the form of the nonlinear differential- 
difference equations which is compatible with these conditions. This proces give rise to 
classes of equations. The list of all nonlinear intcgrable differential-difference equations is 
in principle very long, as any point transformation transforms an integrable equation into 
another one of the same kind. Thus we will consider two equations different only if they 
cannot be related by a point transformation. The class of nonlinear differential-difference 
equations we will study in the following is given by 

Un,t{i) = fn{Un-l{t), Unit), Un+lif)), (1-1) 

where is a complex dependent field expressed in terms of its dependent variables, 

t varying over the complex numbers while n varying over the integers. Eq.(l.l) is a 
differential functional relation which correlates the time evolution of a function calculated 
at the point n to its values in its nearest neighbor points (n + 1, n— 1). A peculiarity of 
the choice of eq.(l.l) is the fact that the right hand side of it is not just a function, i.e. 
it is not the same for all points in the lattice, but for each point of the lattice one has 
an a priori different right hand side. In such a way, for example, by imposing periodicity 
conditions on the dependent field in the lattice variables one is able to rewrite eq.(l.l) 
as a coupled system of nonlinear differential-difference equations. For example if 

/2n = U2n{u2n+l — U2n-l), f2n+l = U2n+2 — U2n, (1-2) 

we can define 

an = U2n, K = U2n+l, (1-3) 

and eq.(l.l) reduces to 

an,t = an{bn-bn-l), bn,t = dn+l - O^n (1-4) 

which by introducing 

~ ^ ) bn — Xn,t 

reduces to the Toda lattice equation 

Xn,tt = e^'^+'-^- - e^"-^-K (1.5) 

Before, in the framework of the formal symmetry approach, only n-independent 
differential-difference equations were considered; the following classes of equations were 
completely classified: 

Un,t = f{Un-l, Un, Un+l) (1.6) 
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(Volterra type equations, see [2]) and 

Un,tt = f{u ) (1-7) 

(Toda type equations, see [3]). 

If in the case of the class of equations (1.6) an equation is defined by a function /, 
in the case of (1.1) we have an infinite set {/„} of a priori quite different functions. So, 
this paper wih be a next step in the classification and in the general theory of the formal 
symmetry approach (readers can find elements of a previous version of the general theory 
only in [3, 5]). 

Section 2 is devoted to the construction of a certain number of conditions necessary 
to prove that an equation of the class (1.1) has generalised symmetries and higher order 
conservation laws. The obtained conditions arc applied in Section 3 to define two main 
classes of nonlinear equations possessing generalized symmetries. The complete classifi- 
cation of this classes of equations is left to a future work. Here one just shows that the 
two classes arc not empty as in one case we have the Toda lattice while in the other the 
discrete version of the Krichever-Novikov equation. 



2. Construction of the classifying conditions 

If equation (1.1) is to represent a nonlinear evolutionary difference equation, then it 
must depend in an essential way from the points (n ± 1), the nearest ncibouring points 
with respect to the point n in which we compute the time evolution. This implies that 
we must add to eq.(l.l) the condition 

dfn -^^0 for any n. (2.1) 



dUn+l dUn- 

Before considering in detail the problem of costructing generalized symmetries to eq.(l.l) 
we will introduce few definitions necessary for the future calculations. 

A function Qn depending on the set of fields Un, for n varying on the lattice, will be 
called a rectricted function and will be denoted by the symbol RF if it is defined on a 
compact support, i.e. if 

), ^ > J, (2.2) 

and i and j are finite integer numbers. 

If there exist, in the range of the possible values of n, values k and m such that 

^^0, ^^7^0, (2.3) 
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then wo say that the function £f„ has a length i — j + 1. For example, Qn could be given 
by the function: 

Qn = nun+i + w„ + [1 + (-l)"]u„_i; 

then i = 1, j = —1 and the length of 5„ is equal to 3 even if only the even functions are 
depending on Un-i- 

Let us define the shift operator D such that 

Dgn{Un+i, Un+j) = gn+l{Un+i+l, Un+j+l)- 

Then we can split the RF into equivalent classes. 
Definition. Two RF 

dn (Wn+ia ) • • • ) '^n+ja ) Sud 6„ [Un-\-if, , ■ ■ ■ , Un-\-jf, ) 

are equivalent 

iff 

an-bn = {D- l)Cn, (2.4) 

where c„ is a RF. 

If, for example, we have Un = Un + Un+i, it is immediate to see that a„ is equivalent to 
a function 6„ = 2u„, as 

an-bn = Un+1 - U„ = {D - 

Let us notice that any function which is equal to a total difference is equivalent to 
zero, i.e. a„ = (D — l)c„ ~ 0. If a RF a„ of length i— j + 1 {i > j) is equivalent to zero, 
then there will exist, by necessity, a RF c„ of the length i — j such that an = {D — l)c„. 
More precisely. 

We see that 

da-n _ dcn+i . 

dUn+i dUn+i 

and consequently 



d^ar, 

0. (2.5) 



dUn+idUn+j 
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In the case i = j, 

O-niUn+i) = C„+l(u„+fci+l, . . . , U„+fc2+l) — Cn{Un+ki, ■ ■ ■ j^n+fe)- 

But as Q^^^ = - Qu^l^^ = for < i, then c„ = rf„(u„+/ci , • • • , Un+i) and consequently 

O-niUn+i) = dn+\{Un+ki+\i ■ ■ ■ ,Un+i+l) — dniUn+ki, ■ ■ ■ ,Un+i)- 

But also g„^"^ = for fci > i so that rf„ cannot depend on Un+k for any k and 
consequently 

= 0. (2.6) 

dUn+i 

So a„ is invariant as, by definitfon, invariant functions are those which depend only on 
n and not on m„. 

We can moreover define the "formal" variational derivative of a RF a„of length 
i — j + 1 as: 

San ^ dak_ ^2 7) 



This quantity is strictly related to the notion of variational derivative, and this is the 

San 
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reason for its name. If a„ is linear, then |^ is an invariant function, but if it is nonlinear 



then ^ 

-T-^ = 9n{Un+N , Uu-n), (2.8) 

where for some k ^ and for some m „ ^"^ ^ 0. It is immediate to prove 

that if a„ is a RF which is equivalent to zero, i.e. it is total difference, i.e. it is the 
difference of an RF calculated in two neighboring points of the lattice, then the formal 
variational derivative of a„ is zero. One can also prove that if |^ = 0, then a„ is 
equivalent to zero. In fact, using eq.(2.7) we have 



dUn+idUn+j 

which implies that 

0'n = bn{Un+i, U„+j+l) + C„(u„+i_i, Un+j) ^ dn{Un+i-l, Un+j), 

i.e. a„ is equivalent to a RF of length i — j. Carrying out recursively this reasoning, we 
arrive at the conclusion that 

a„ ~ Fn{un) with = 0, 
5 



i.e. a„ must bo equivalent to an n-dependent constant (or invariant function) which is 
always equivalent to zero. 

Given a nonlinear chain (1.1) we will say that the RF gn{un+i, Un+j) is a 

generalized (or higher) local symmetry of our equation of order i (more precisely, of left 
order i) iff 

Un,T = 9n{Un+i, Un+j) (2.9a) 

is compatible with (1.1), i.e. iff 

dtdr{Un)=drdt{Un). (2.96) 

Explicitating this condition we get: 

Ctdn = OtJu = t; Un+1,T + 7^ Un,T + ^ Wn-l,r 

OUn+1 OUn OUn-1 



dfn _|_ dfn ^ dfn 



dUn+l dUn dUn-1 



9n — fn9n, 



I.e. 

{dt - f:)9n = 0, (2.10) 
where by /* we mean the Prechet derivative of the function /„, which is given by: 

f* = ^^-^ n + ^ + ^-^^ _n-i = /■(!) n + fio) . /•(-!) n-i (o u) 

Eq.(2.10) is an equation for £f„ once the function /„ is given, an equation for the sym- 
metries. In this work we limit ourselves to local symmetries, i.e. symmetries which are 
RF. A nonlocal extention can be carried out by introducing, for example, a new field Vn'- 
{D - l)v„ = Un, i.e. Vn = -Yl'jLn'^j '^^ '^n = (compare [5]). Extension in 

such a direction will be carried out in future work. 

Given a symmetry we can construct a new symmetry by applying a recursive oper- 
ator, i.e. an operator which transforms a symmetry into another one. An operator 

m 

Ln= ^n\t)Di (2.12) 



will be a recursive operator for eq.(l.l) when 

9n = Ln9n (2.13) 
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is a new generalized symmetry associated to (1.1). This condition implies that 

^(L„)=L„,t-[/*,i„]=0. (2.14) 



From (2.10) it follows that 

A{g*J=dr{f;). (2.15) 

In fact, defining 

Bn = dt{gn) = Q ^" /w+fc) 

^ OUn+k 

we have 

= a /^"a « Jn+kJ-> + / ^-R -7. ^ 

m m,k m,k 

= dt{g:) + g*j:, (2.16) 

and thus cq.(2.15) is obtained by introducing (2.16) into the Frechet derivative of (2.10). 

Eq.(2.15) implies that, as the rhs is an operator of the order 1 (see (2.11)). the 
highest terms in the Ihs must be zero. Thus we can define an approximate symmetry of 
order i and length m, i.e. an operator 

k—i—m+l 

such that the highest m terms of the operator 

i+l 



A{Gn)= Yl ^n^D' 



k=i—m 

are zeroes. Taking into account eq.(2.15), we find that we must have i — m + 2 > 1 if the 
equation 

A{Gn) = (2.17) 

is to be true and G„ is to be an approximation to the Frechet derivative of gn, the local 

generalized symmetry of our equation. 

From this results we can derive the first integrability condition which can be stated 
in the following Theorem: 



Theorem 1. If eq.(l.l) has a generalized local symmetry of order i > 2, then it must 
have a conservation law given by 

ailog/W = (D-l)gW, (2.18) 

where is an RF. 

The proof of this Theorem is contained in Appendix A. 

In this way we have showed the existence of the first canonical conservation law. 
The next canonical conservation laws could be obtained in the same way, by assuming 
a higher symmetry and thus being able to consider an approximate symmetry of higher 
length. These canonical conservation laws would, however, be very complicated (they 
will depend on the order of the generalised symmetry) and very difficult to reduce to 
simple expressions (not depending on this order). So we prefer to follow an alternative 
way which requires the existence of two higher symmetries. 

This procedure can be carried out, as we already know one canonical conservation 
law. Let us now assume that there exist two RF gn and ^„ which generate symmetries 
of order i and i + 1 respectively. In correspondence to the symmetries, wc can construct 
the approximate symmetries Gn and Gn of orders i and i + 1 respectively, and from (2.1) 
Qn^ and gn~^^^ will be different from zero for all n (see (A3), Appendix A). Starting from 
Gn and G„, we can construct the operator 

Gn = G-^Gn. (2.19) 

As from (2.14) we have 

A{G-^) = -G-^A{Gn)G-\ A{LnKn) = A{Ln)Kn + LnA{Kn), 

we obtain 

A{Gn) = G-'[-A{Gn)Gn + A(G„)]. (2.20) 

Let's notice that, as Gn is an approximate symmetry, its inverse will be an operator char- 
acterized by an infinite sum. Consequently G„ , thought being an approximate symmetry 
of order 1 and length i (the lowest of the two lengths of G„ and Gn) is represented by an 
infinite sum. This shows that under the hypothesis that two local higher simmetrics exist, 
we can reduce ourselves to consider approximate symmetries of order 1. In such a way 
Qn^ = fn^ and for the following simple formula can be obtained: qn^ = gn^ — fn\ 
We can now state the following Theorem proved in Appendix B: 

Theorem 2. If the equation (1.1) with the conditions (2.1) has two generalized local 
symmetries of order i and i + with i > 4, then the following canonical conservation 
laws must be true: 

atpW = (Z?-l)gW (fc = 1,2,3), 
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(1) 1 ^/n (2) (1) , ^fn 



Pn -Qn +^{Pn ) + Qu^^^ dUn ' ^ ^ ^ 



where q^'^ (fc = 1,2,3) are some RFs. 

As we have seen up to now, if eq.(l.l) has local generalised symmetries of high 
enough order, we can construct a few canonical conservation laws depending on the 
function at the rhs of (1.1). 

One can divide the conservation laws into conjugacy classes under an equivalence 
condition. Two conservation laws 

Pn,t = {D - l)g„, r„_t = (-D - l)s„ 

are said to be equivalent if 

Pn ~ r„. (2.22) 

Then we can say that the local conservation law is trivial it Pn ~ 0. it Pn ~ fn {un), with 
^ at least for some n, then we say that the conservation law is of zeroth order while 

if 

Pn ~ fn {Un+N, • • • , Un), N > 0, 

and 

dUndUn+N 

for at least some n, wo say that the conservation law is of order N. 

An alternative way to define equivalence classes of local conservation laws is via the 
formal variational derivative. If the local conservation law is trivial, then |^ = 0, if it is 

of zeroth order, then |^ = ^ni'^^n) 7^ for at least some n while if it is of order N, then 
Spn . , . 

=Pn{Un+N, Un, Un-N), 

OUn 

where 

dUn+N ' dUn-N 

for at least some n. 

Let us denote by p„ the formal variational derivative of the density p„ of a local 
conservation law, i.e. 

Pn = (2.23) 

SUn 
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For any conserved density p„ we derive 

Pn,t ~ Pnfn ~ (2.24) 

by direct calculation; then, by carrying out the formal variational derivative of eq.(2.24), 
taking into account that in a summation 

dpn+k dpn 



dUn dUn+k 

we get that the formal variational derivative p„ of a conserved density p„ satisfies the 
following equation: 

(5t + /DPn = 0, (2.25) 
where the transposed Prechet derivative of /„ is given by 

Let us consider the Prechet derivative of p„ for a local conservation law of order N. 
In such a case, we have 

Pi = E P'^'D'^ P'n' = (2-27) 

k=-N 

Let us construct the following operator: 

B{p*n) = P*n,t +P:f: + ffpl = E (2-28) 

k 

where are some RFs. From eq.(2.28) we get that 

= P^} + TMi'ftP + f!^^i^f)- (2-29) 

j 

By differentiating eq.(2.25) with respect to Un+k, we can rewrite eq.(2.29), after a long 
but straightforward calculation, in the form 

10 



and thus prove that, as /„ depends just on u„ and u„±i, bn are different from zero only 
for -2 < fc < 2. 

In such a way, for a sufficiently high order conserved density p„, we can require that 

Bipi) = (2.31) 

is approximately solved. If the first t77. < iV — 1 terms of the Frechet derivative of Pn 
satisiy eq.(2.31), then we say that we have an approximate conserved density of order N 
and length m. 

Taking all the results up to now obtained into account, we can state the following 
Theorem which will be proved in Appendix C: 

Theorem 3. If the chain (1.1) with conditions (2.1) has a conservation law of order 
iV > 3, and the condition (2.18) is satisfied, then the following relations must take place: 

r«=(Z)-l)s« (fc=l,2) (2.32a) 

with 

ril' = log[-fi'^/fi-% r(^) = s^:l + 2/(°), (2.326) 

where Sn^ are RFs. 

We will end this Section by adding a few comments. First of all we can notice that the 
request of existence of nontrivial local conservation laws is more restrictive than that of 
existence of symmetries. In fact, there are many instances in which there exist generalized 
symmetries, but not nontrivial conservation laws. This is the case, for example, of the 
c-integrable equations, i.e. such nonlinear equations which are transformed into a linear 
equation. 

If one compares Theorem 1 and Theorem 2, one can think that among conditions 
(2.18) and (2.21) with A; = 2,3 there is a difference of importance, as conditions (2.21) 
require the existence of two generalized symmetries while for condition (2.18) only one 
generalized symmetry is sufficient. However, we could obtain conditions (2.21) with 
fc = 2, 3, assuming that only one symmetry of an order i > A exists, but calculations in 
the proof would be much more difficult. 

All cases cosidered here required that just fn^ ^ when we derived conditions 
(2.21). An analogous set of conditions could be derived if we request that just fr-i~^^ ^ 
for all n. They can be derived in a straitforward way considering expansions in negative 
powers of D, instead of positive as we have done here. This derivation is left to the readers 
as an exercise. This set of conditions also will have the form of canonical conservation 
laws: 

dtpi'^ = {D-l)e\ (2.33) 
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and conserved densities will be symmetric to ones of (2.21). For example 



p«=log-^. (2.34) 

OUn-l 

Let us notice moreover that if Hn^ and Hn^ are two solutions of (2.31) of different 
order, the operator 

K„ = (H^'^r'Hi^^ (2.35) 

satisfies (2.14) and thus it is a recursive operator. Consequently, if we start from two 
approximate solutions of (2.31), i.e. two Frechet derivatives of formal variational deriva- 
tives of conserved densities, wc can, using (2.31), get an approximate symmetry. So, one 
can derive all the conditions (2.21), (2.32), (2.33), assuming the existence of two higher 
order local conservation laws. In the case of conditions (2.33), one should use the same 
formula (2.35), but {Hn^)~^ will be a series in positive powers of the shift operator D: 

+ 00 



k=-N 

This proves the statement written before that conservation laws are "more fundamental" 
than symmetries. 

If we compare conditions (2.21), (2.32), (2.33), we can see that, for example, 

rW=log(-l)+pW-pW, 

i.e. the first of conditions (2.32) means that pn ^ ~ pn \ i.e. first canonical conservation 
laws of (2.21) and (2.33) are equivalent. The same result could be obtained for the second 
condition of (2.32). In particular, the set of conditions (2.33) can be derived starting 
from conditions (2.21) and (2.32). However, these conditions are of great importance in 
themselves, as there might be equations of interest which satisfy (2.21) and (2.33), but 
not (2.32). 

The solution of the conditions, be those obtained by requesting existence of the 
generalised symmetries or those of local conservation laws, provide the highest order 
coefficients of the Frechet derivative of a symmetry or of the formal variational derivative 
of a conserved density. Those coefficients are the building blocks for the reconstruction of 
the symmetries or of the formal variational derivatives of the conserved densities. In fact. 



limiting ourselves to the case of symmetries for the sake of concreteness, gi^^ = 



for k = i,i — 1,.... Thus the knowledge g„ ' for a few values of k gives a set of partial 
differential equations for gn in terms of its variable whose solution provides the needed 
symmetry. The same way allows us to reconstruct variational derivatives of conserved 
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densities. There is, however, a more direct way to obtain conserved densities. In fact, 
if we know highest coefRcients of the first order solution L„ of eq.(2.14), we can obtain 
several conserved densities by the following formula: 

p«=res(L^J (j = l, 2, ,,,) 

(see Appendix B). Moreover, the solution of the conditions (2.21) leads to three local 
conservation laws in explicit form. 

3. Introduction to chain classification 

In this Section, using just two conditions, we give two Theorems which split eq,(l.l) 
into two classes of nonlinear equations depending on a limited number of variables. The 
complete classification using all other conditions will be presented in the second part 
of this work. Here, however, we present at the end examples of nonlinear evolution 
equations on the lattice belonging to each of the obtained classes. 

3.1 Main classification Theorems 

Let us apply the first of the conditions (2.32) to eq. (1.1), (2.1), i.e. 

r(i) = log[-/(^V4-'^] = {D- l)s«. (3.1) 

In such a case, 



^n^ = ?i^Hw.n+l,Wn,W.n-l) 

I 

Introducing the function 



and consequently Sn^ = Sn\un,Un-i) is defined uniquely if Sn^ is to be a RF. 



Vn = yn{Un+l,Un) = exps^^|i, 

we get from (3.1) that y„ must satisfy the following difference equation: 

ynfi-'^+yn-ifi'^ = o (3.2) 

where the function y„ is, by definition, in a generic point of the space of its independent 
variables, always different from zero. Eq.(3.2) is a constraining equation for /„ which 
thus provide a reduction of /„ in terms of new symmetry variables. In fact, defining a 
function t;„(u„+i, m„_i) so that 

9vn dv„ 
Vn = ) Vn-i = , (3.3) 
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we find that (3.2) is identically satisfied if /„(u„+i, u„, u„_i) is written as 

fn = 'Pn{Un,Vn)- (3.4) 



Then, as 



we get that 



Vn. (3.5) 
Eq.(3.3) can be "solved" by introducing a function z„ = 2„(u„_|_i, u„) such that 

Vn = . I (3.6) 

and 

d 

Vn = ^^(-2" ~ Zn-l)- (3.7) 

We can now state the following Theorem: 

Theorem 1. The condition (3.1) is satisfied iff the rhs of (1.1), taking into account 
conditions (2.1), can be represented as a function (3.4) with w„ given by (3.7) in term of 
a function ^„(u„+i, u„) such that 

-^^^Q Vn. (3.8) 

The function ipn will depend in an essential way from the variable Vn (see (3.5)). 

The proof is already contained in the results presented before; (3.8) is given by (3.6), 
as Dn is always different from zero. In such a way, we have reduced the problem of finding 
a set of functions of three variables to that of finding two sets of functions of two variables 
only. Even if w„ depends on u„, the variables u„ and u„ may be taken to be independent 
variables, as Un is always different from zero. 

We can now pass to consider the first of the canonical conservation laws given by 
(2.21), i.e. 

Pi^]~0, p(^)=log/W=log(^y„). (3.9) 
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Let us explicitate (3.9); we have 

= fn+1 + Jn + r. In ~ Pn t ~ 0' (3-10) 

dUn+\ OUn OUn 

As is a RF equivalent to zero, we have 



dUn+2dUn-l dUn+ldUn-1 dUn+2 dUn+2dUn dUn-l 

Prom the definition of pl^^ in terms of and y„, we get 



Vn. (3.11) 



V dvn+ij dv„+i dvl \ dvnj dv„ 

To obtain (3.12), we have used the fact that i/n and are both different from zero for 
any n. Eq.(3.12) implies that the function <fin{un,Vn) must satisfy the following ordinary 
differential equation: 



diPn 2 



dv„ 



eip^ + an{un)'Pn + bn{u„), (3.13) 



whore e is an arbitrary n-independent constant, and a„ and 6„ are two arbitrary in- 
dependent sets of functions. 

Eq.(3.13) implies that wc have two main different cases, characterised by the form 
of the dependence of (pn on i>„: the case £ ^ and the case £ = 0. If £ ^ and, for 
example, 

A„(u„) =al- 4ebn Vn, (3.14) 



then 



A^^ {aI/^ \ an 

n = tanh I —^{Kn{Un) -Vn)\ - ^, (3.15) 

where Kn is an arbitrary integration M„-dependent set of functions. If 

£7^0, A„=0 Vn, 

we are led to 

s{Kn{Un)-Vn) 2e 
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In general, if we consider the case e 7^ 0, the function (pn may have (at least theoretically) 
the form (3.15) for some n and the form (3.16) for some other n. At the end of this Section 
we will write down an integrable example of this kind. 

Let us consider the case £ = 0. Eq.(3.13) is equivalent to the following ordinary 
differential equation: 

= exp[a„('u„)?;„ + c„(w„)], (3-17) 

where c„ is an arbitrary u„-dependent set of functions. It is easy to show that the 
function (p„ must have the form: 

= an{Un) exp(a„ti„) + (3n{Un)Vn + ln{Un). (3.18) 

Using (3.13) or the equivalent equation (3.17), we can obtain conditions for coefficients 
ocni Pn, In- In particular, a„/3„ = for any n. Below readers can see an example such 
that 

where P„ and Qn are two orthogonal projection operators defined by: 

Pn = Qn = \ ■ (3-19) 

It is clear that the equation Q!„/?„ = has not only such the solution; the Toda and 
Volterra equations correspond to the solution: q;„ = 0, /3fi 7^ for any n. 

Eqs. (3.15), (3.16), (3.18) provide the two classes of nonlinear differential-difference 
equations of type (1.1) which satisfy the simplest integrability condition of (2.32) and 
partially the condition (2.21) with A; = 1. The other conditions (2.21) and (2.32) with 
respect to an arbitrary point transformation will reduce almost completely the classes of 
equations to a few equations only. 



3.2 Examples 

Let us consider an example of equation belonging to the first class, more precisely 
to the class 

e^O, An^O Vn. 

In such a case we have to fix the arbitrary functions iir„(u„), a„(u„) and 6„(u„), the 
arbitrary constant e, and to define the function Vn as function of Un, Un-i to get 

a well defined nonlinear evolution chain equation. 
Let us set: 

Kn = 0, a„ = 0, hn = ^\p{un)r{un) - q^{un)], ^ = ^' (^-20) 
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where 



p(u„) = q;m^ + 2/3u„ + 7, (3.21a) 
r{un) = 7?4 + 25u„ + w, (3.216) 
q{un) = Pul + XUn + 6. (3.21c) 



Then we can define 



2 , , f p{Un)Un+lUn-l+q{Un){Un+i+Un-l)+r{u„)\ 

Vn = — —7= arctann I r^=^- ] (3.22) 

V-26„ V V-26„ {Un+l - Un-l) J 



which corresponds to define 

Un 



• + 2q{Un)Un+l + r{Un) ' 

Inserting (3.20), (3.22) into (3.15), we get the following evolutionary nonlinear chain: 

P{Un)Un+lUn-l + q{Un)(Un+l + Un-l) + r{Un) 

Un,t = • {6.Z6) 

Un+l — Un-l 

Eq.(3.23) depends on 6 arbitrary complex constants (sec (3.21)) and is invariant 
under linear-fractional transformations (only coefficients of the polinomials p, q, r are 
changed). It was obtained for the first time in [2] when classifying discrete evolutionary 
equations of the form: 

Un,t = fiUn+1, Un, Un-l). (3.24) 

It satisfies all the five integrability conditions, has an infinite set of higher local conserva- 
tion laws and should have an infinite set of generalized symmetries (but nobody has yet 
proved it). It is the only example of nonlinear chain of the form (3.24) which cannot be 
reduced to the Toda or Volterra equations by Miura transformations. By carrying out 
the continuous limit, as one does for the Volterra equation to obtain the Korteweg-de 
Vries equation, we get the well-known Krichever-Novikov equation [8] : 

Wt = u..^^-|— + ^^, (3.25) 

Z Ux Ux 

where R{u) is an arbitrary 4th degree polinomial of its argument with constant coeffi- 
cients. 

In the second case £ = 0, let us restrict ourselves to the subcase a„ = for any n 
(see (3.18)) and consider functions of the form: 

= Pn{Un)Vn + IniUn)- (3.26) 
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Let 

7„ = 0, = Un+l - Un-1, (3.27) 

but j3n be an arbitrary nonzero (for all n) function of its argument. In such a case, 

fn = f3n{Un){Un+l - Un-l)- (3.28) 

The first canonical conservation law implies that 

dtlogPn ~ 0, 

i.e. 

f3n(Un+l - Un-l) ~ /3>n+l " P'n+lUn ~ (3.29) 

from which we derive that (3!^ must be an n-independent constant: 

= Aul + + C„. (3.30) 
Inserting this result into (3.29), we get that 

(B„_i - B„+i)u„ ~ 0, i.e. Bn=B + (-1)"B. 

Consequently, the first canonical conservation law reads: 

where 

q^^^ = 2AUnUn-l + B{Un + Un-l) - B[{-l)"Un + (- 1)"-1m„-i] . (3.31) 

Introducing (3.28) and (3.31) into the second canonical conservation law, we get, 
after a straightforward but lengthty calculation, that 

pi'l - {Cn-l " Cn+l){Aul + BnUn) ~ 0. 

We can see, in particular, that ii A^Q, then 

(7„ = C+(-l)"C. (3.32) 

Using the last canonical conservation law, one can prove that C„ always must have the 
form (3.32). So, we are led to the chain 

Un,t = [Aul + {B+ (-l)"^)u„ + C+ (-l)"(7](u„+i - Un-l) (3.33) 
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with 5 arbitrary complex constants. 

By direct calculations, one can then prove that this nonlinear equation (3.33), for 
any choice of the arbitrary constants, satisfies all three canonical conservation laws (2.21) 
and both additional conditions (2.32). By obvious point transformations, we can reduce 
any nonlinear chain of the form (3.33) to one of the following chains: 

Un,t = 0n{Un){Un+l - Un-l), (3.34) 

where 



(the Toda chain), 

(the Volterra equation) or 

and in the last case, 



Pn = PnUn + Qn (3.35) 
(in = Un (3.36) 
Pn = Cn- ul, (3.37) 



Cn = l, Cn = or Cn = Pn (3.38) 

(three modifications of the Volterra equation). Here P„ and Qn are the orthogonal pro- 
jection operators (3.19). Unlike the discrete version of the Krichever-Novikov equation 
(3.23), (3.21), all the chains (3.33) can be reduced to the Toda chain by Miura transfor- 
mations. For example, in the case of the Volterra equation, we have the transformation 

Un = PnUn+lUn + QniUn+l+Un). (3.39) 

Formula (3.39) brings any solution m„ of the Volterra equation into a solution £(„ of the 
Toda chain. Transformations of the modified Volterra equations (3.34), (3.37), (3.38) 
into the Volterra equation are given by the formula 

Un = {Cn + Un){Cn+l - Wn+l) (3.40) 

in all the three cases. Transformations (3.39), (3.40) together with point transformations 
enable one to prove that any nonlinear chain of the form (3.33) possesses local conser- 
vation laws of an arbitrary high order. This means, in particular, that the chains (3.33) 
satisfy not only the classifying conditions (2.21), (2.32) but also all conditions we could 
derive using approximate symmetries and conserved densities. 
The last class of equations we discuss here is given by: 

Un,t = Fn{Un+l-Un-l), F^^ = AF^ + B^Fn + Cn, (3.41) 
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where 



B„ = B + (-l)"B, C„ = C + (-l)"C, 



and A, B, B, C, C are arbitrary complex constants. These equations have infinite sets of 
local conservation laws and satisfy all the five classifying conditions because are reduced 
to chains of the form (3.33) by the following transformation: 

Un = Fn{Un+l - U„_i). 

In the case 

Fn{z) = Pntanh^; + QnZ'^ (3.42) 
one has the differential equation 

F' = P — F"^ 

i.e. function (3.42) defines an equation of the form (3.41) corresponding to the last of 
the modified Voltcrra equations (sec (3.34), (3.37), (3.38)). So, wo have obtained an 
example of the chain with the rhs which can be represented as function (y9„(u„, Vn) (here 
Vn = Un+i — Un-i) being tanhv„ for odd n and for even n. The second example of 
this kind is given by: 

Fn{z) = P„expz + QnZ. 

In this case, 

i.e. the equation of the form (3.41) is reduced to the Toda chain. 
Appendix A: Proof of Theorem 1 

For a sufficiently high order symmetry, i.e. i >> 1, the highest terms of 

will satisfy the following equation: 

i:9>' + E E 9it^f'::iD'^- -EE fi"''9'nlmD'^"' = 0' (^1) 

1=2 k m= — 1 k m=— 1 

where the sum is over those k such that fc + m > 1, as otherwise the Ihs of (Al) is no 
more zero. In (Al) the coefficients of any power of D must vanish; so the coefficients of 

reads: 

^^^Vl^i = f^^^9%. {A2) 
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As, due to (2.1), fn^ ^ Vn, we have 



n ft'' 



k=n 



where we have, with no rectriction, set to unit the "integration" constant. 

Let us consider now the coefficient of this conies from more than one term (fc = i, 
TO = Oorfc = i — l,m = l) and involves the time evolution of Qn^ . It can be cast in the 
following form: 



(i) 



(»-l) 
9n+l 



n+i — 1 



(i-i) 

_ en - n V 

(1) „n+i-2 „(1) 
llfe=n+l/fe llfe=n Jk 



(0) 



k=n 



from which we derive: 



(i-l) n+i-1 

9n .(0) 

„(1) Jk 
_llfc=n Jk k=n 



Introducing (A3) into the Ihs of (A5), we can write it as: 



n-\-i—l n-\-i—l 



9tiog n fk = E dtlogfl'> ^idtlogfi'^ ^0 c.v.d. 



k=n 



k=n 



{AA) 



(A5) 



Appendix B: Proof ot Theorem 2 

As a consequence of Theorem 1, we can assert that we have an approximate sym- 
metry of order i = 1 



where from (A3) 
Instead of (A5) we have 



G„ = g^)D + + gi-'^D-^ + gi-^)D-^ + ...., 



a,log/« = (Z)-l)(^)-/(°)) 



{Bl) 

m 
m 



(see the coefficient of D in the equation A{Gn) = with A defined by (2.14)). Conse- 
quently, the function of the rhs of the first canonical conservation law is: 
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from which we derive 

Let us now consider the coefficient of -D° in the equation A{Gn) = 0: 

dtgl^^ = {D- l)[£lAgk-'^ - ft'^)]- {B5) 
Eq.(B5) is the second canonical conservation law with the Ihs given by (B4): 

Prom (B6) we get 

To get the last relation, it is simpler to use the following Lemma: 

Lemma. If = hn^D^ + hn~^^D^~^ + ■ ■ ■ is an approximate symmetry, which 
satisfies the first m> i + 2 terms of the equation Hn,t = Hn], then 

res(if„) = /i^o) (B8) 

will be a conserved density. 
In fact, 

dtresiHn) = rcs(ff„,0 = TCs[f*,Hn]. {B9) 
The coefficient of D° of [/*, Hn] can be obtained only from terms of the type 

which are equivalent to zero. 

As any power of an approximate symmetry is also an approximate symmetry of the 
same lengh, we can construct a new conserved density, calculating the residue of G^. In 
such a case, after a long but straightforward calculation, we get 

resid) = yesiigi^^D + + gi-'^D-^ + ■■■)'] = 

As, from the previous Lemma, if dtresC^ ~ 0, so will be dtPn^ ~ c.v.d. 
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Appendix C: Proof of Theorem 3 

Let us assume that we have a solution of 

dtPn + Pnfn + ffPn = 0, (CI) 

where Pn is an approximate conserved density 

JV 

Pn= E Pn^^' ^^2) 
fe=Ar-m+l 

of order iV > 3 and length m > 2. In such a case, introducing (C2) into (CI), the 
coefficient of D^^^ reads: 

As /n^^' for any n, we can prove that p^^^ ^ for any n. Then we get 

and thus, by taking the logarithm of both sides, we are led to 

r«=log[-/W//(-i)]. (C5) 

From (C4) and (C5) we get 
where s„ is such that 

+ fk-^^Sn+i = 0. (C7) 
The coefiicient at order of (CI) gives: 

d, l0g(s„+i/«i • • • f^nU-l) + f^n^ + f^aU ~ 

from which 

(AT - 1) dt log /(I) + a* log s„ + 2/(0) ^ 0. 

As condition (2.18) is satisfied, and logs„ = Sn^ (compare (C7) and the first of the 
conditions (2.32)), we are led to 

r(^)=a,a«+2/10) cv.d. 
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